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Performance studies of the two-step multiboson algorithm in compact 
lattice QED * 

I.L. Bogolubsky'", V.K. Mitrjushkin^, I. Montvay'^, M. Muller-Preussker'= and N.V. Zverev" 
^Joint Institute for Nuclear Research, 141980 Dubna, Russia 
''Deutsches Elektronen-Synchrotron, 22603 Hamburg, Germany 
"^Humboldt-Universitat zu Berlin, Institut fiir Physik, 10115 Berlin, Germany 

The performance of the two-step multiboson (TSMB) algorithm is investigated in comparison with the hybrid 
Monte Carlo (HMC) method for compact lattice QED with standard Wilson fermions both in the Coulomb 
and confinement phases. The restriction to QED allows for extensive measurements of autocorrelation times. 
Preliminary results show that the TSMB algorithm is at least competitive with standard HMC. 



Compact lattice 4D QED with Nf = 2 dynami- 
cal Wilson fermions has been numerically studied 
jl] by means of the HMC algorithm |Q . However, 
the use of this algorithm within the confinement 
phase in the chiral limit k — > Kc{P) as well as 
in the strong coupling phase with k > Kc (the 
so-called 'Aoki phase') becomes cumbersome [|j 
because of large average condition numbers 
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Here (Amax) and (Amin), respectively, are the 
largest and smallest average eigenvalues of the 
even-odd preconditioned Wilson fermion matrix 
A4^A4. Another open principal problem for the 
HMC method is the simulation of an odd number 
of fermion flavours. Hence, an alternative to the 
HMC method is desired. The same arguments 
hold also for lattice QCD and related models. 

An alternative is provided by the hermitian 
two-step multiboson (TSMB) algorithm invented 
recently by one of the present authors (I.M.) 
m. It consists of the following basic ingredients. 
First, a polynomial approximation of the inverse 
fermion determinant is introduced. Second, 
in order to reduce the computational costs [||-^ 
a noisy correction accept-reject step is employed. 
At the end one performs an auxiliary reweighting 
step for correcting the measured observables. 

The TSMB algorithm has been successfully ap- 
plied to a supersymmetric model 0]. Here we 



would like to investigate the efhciency of this al- 
gorithm in comparison with the standard HMC 
one. For this aim compact lattice QED with Wil- 
son fermions provides a valuable test ground, be- 
cause its fermionic properties - in particular in 
the confinement phase close to the chiral limit - 
resemble very much those of lattice QCD with 
Wilson fermions In this paper we consider 
always the case Nj ~ 2. 

First of all, we want to establish a strategy for 
appropriately choosing the technical parameters 
for both the HMC and the TSMB algorithms. In 
order to achieve a reasonable acceptance rate for 
HMC we adjust the number of time steps Nr and 
the time step size At in the Hamilton dynamics 
according to the rule 
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where C is defined in (P, and V = x N4 de- 
notes the lattice volume. The stopping criteria 
for the conjugate gradient (CG) invertor in the 
Hamilton dynamics i5hd and for the accept-reject 
step Sacc are fixed in accordance with pd| 
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As a consequence the average number of CG it- 
erations can be estimated to 
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*Talk presented by N.V. Zverev at LATTICE 2001, Berlin. 



The TSMB algorithm uses polynomials within 
the least squares integral approximation of the 
function x'^^l"^ over the interval [e. A] (see 
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|T^). The first (crude) polynomial of order ni 
is used in the first multiboson step. The second 
(correcting) polynomial of order n2 is required for 
the noisy correction step. The third polynomial of 
order 77,3 approximates the inverse square root of 
the second polynomial. And the fourth (finest) 
polynomial of the order 714 is necessary for the 
reweighting step. These parameters are chosen 
according to a prescription in ]l3t : 

77ioc(VC)'/'lnl^, 772 ^ (A^^^^)), (5) 

713 ~ (1.2 1.4)772, 714 > 772, (6) 

A~ (1.2— 1.4)(A„,ax), e~0.5(An.in). (7) 

For 7ii we have taken into account the empiri- 
cal observation ||l^ that the dependence on the 
condition number becomes effectively weaker for 
the least squares first polynomial compared with 
polynomials in the Chebyshev approximation re- 
sulting in a replacement C — \/C- 

Now let us come to the theoretical cost analysis. 
We define the performance of an algorithm by 

P = NopcrTint, (8) 

where Nopcr is the total number of operations 
per update. Tint denotes the integrated autocor- 
relation time for a given observable. The num- 
bers of operations in the HMC and in the TSMB 
(712 2> 771 ) algorithms are given as follows 

N™C^F(7V(C«))7V., N™ocF772. (9) 

For simplicity instead of the value Nopcr we used 
the CPU time required for one sweep tcpv- The 
integrated autocorrelation times for the HMC and 
the TSMB methods can be estimated according 
to arguments given in |lO| and [^,^ 

^HMC ^ (i\T A^^-2 ^TSMB ^ ^ / /7\l/2 /-.r.^ 
Vt OC (A^At) , Tj^j (X7li(VC)' • (10) 

Here the effective dependence on the condition 
number ^ (|l|) for the least squares first polynomial 
has been taken into account, too. Thus, the ex- 
pected theoretical gain G of the TSMB algorithm 
over the HMC method is evaluated according to 
i), (i, (I), (M) as follows 

pHMC yl/i 
G=pTSMB«7^- (11) 



phase 


Nr At (A^^^'^)) 


Coulomb 
confinement 


40 0.025 36.0(2) 
10 0.01 500(2) 



Table 1 

HMC parameters for Coulomb and confinement 
phases. Stopping criteria are Jhd = 10^"^ and 





ni 


"2 


"3 


774 


e 


A 


(a) 


6 


36 


48 


200 


0.025 


2.5 


(b) 


50 


360 


450 


500 


0.000225 
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Table 2 

TSMB parameters for Coulomb (a) and confine- 
ment (b) phases. 



Now we turn to the numerical investigation of 
the performance of HMC and TSMB algorithms 
both in the Coulomb and confinement phases. 
The lattice size is 6'^ x 12 and time-antiperiodic 
boundary conditions for fermions are employed. 
In the Coulomb phase we have chosen the point 
(/3, 7t) = (2.0,0.130), whereas in the confinement 
phase (/3, n) = (0.0, 0.238) quite close to the chi- 
ral limit. For the HMC method the number of 
time steps A^r and the time step size At selected 
according to the prescription (^) are presented in 
Table |l|. The stopping criteria for the CG method 
according to (^) provide an average number of 
CG inversion steps (in the Hamiltonian dynam- 
ics) {N^'~^'^'>) as shown in the same Table. In case 
of the TSMB algorithm the polynomial parame- 
ters are chosen according to the recipe (|l-0) and 
presented in Table ||. The knowledge of these 
parameters requires information about the small- 
est (Ainin) and largest (A,„ax) average eigenvalues 
of the even-odd preconditioned fermion matrix 
Ai^M, tunable to the dynamical fermion case 
from the quenched one (see Table I). 

The performance of the dynamical fermion al- 
gorithms was studied with respect to the aver- 
age values of the following gauge invariant ob- 
servables O: the mean gauge plaquette energy 
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quenched 


dynamical 




(•^min) (-^max) 


(•^min) (-^max) 


(a) 
(b) 


.065(1) 1.60(1) 
.0010(1) 6.78(1) 


.13(1) 1.63(1) 
.0005(1) 6.59(1) 



Table 3 

Minimal (Amin) and maximal (Amax) average 
eigenvalue of the even-odd preconditioned Wil- 
son fermion matrix A4^A4 for both Coulomb (a) 
and confinement (b) phases in quenched and dy- 
namical cases, the latter determined from HMC. 





(Eg) 




(n) 




Coulomb 




0.1332(1) 


0.9381(1) 


1.378(1) 




0.1331(1) 


0.9379(1) 


1.376(1) 


_HMC 
Mnt 


3.2(3) 


2.0(2) 


25(4) 


_TSMB 
'int 


3.0(3) 


2.8(2) 


50(8) 


G 


1.7(2) 


1.2(2) 


0.8(2) 




Confinement 




0.939(1) 


0.95(1) 


13.9(2) 




0.938(1) 


0.96(1) 


13.7(2) 


_HMC 
'int 


65(7) 


60(7) 


35(5) 


_TSMB 
'int 


120(20) 


125(15) 


45(5) 


G 


0.5(1) 


0.5(1) 


0.7(1) 



Table 4 

Performance of HMC and TSMB algorithms in 
the Coulomb and confinement phases. 



(Ea) = (1 — Ret/p) the scalar condensate (i/'^A) 
and the pion norm (11) = (('075'0)^)- We mea- 
sured the corresponding integrated autocorrela- 
tion times Tint and also the gain G (^ of the 
TSMB over the HMC method. 

The main result of our investigation is pre- 
sented in Table ^. The statistics in our case 
(0(10'*) measurements) was sufficient to evalu- 
ate the integrated autocorrelation time. The 
reweighting procedure in the TSMB case has not 
changed significantly the results obtained using 
only the main two steps. We may conclude that 
the performance of the two algorithms is compa- 
rable for the chosen parameters. Further opti- 



mization of TSMB is possible by increasing the 
number of local gauge updates from one (in our 
case) to a larger number proportional to r^f^^ 



in (10). Moreover, the polynomial degrees ni and 
712 could be taken smaller leaving more space for 
error correction by reweighting. This shows that 
for an even number of fermion flavours the TSMB 
algorithm is at least competitive with HMC. 

This research has been supported by the grad- 
uate college DFG-GK 271, by the RFRB grant 
99-01-01230 and by the JINR Heisenberg-Landau 
program. 
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